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ABSTRACT 


Much critical equipment at Kennedy Space Center is plumbing for the transport of liquid Hydrogen or 
liquid Oxygen. Every piece of hardware for handling such a hazardous cryogen is subject to testing prior to 
installation and use. 

Safe, realistic testing of all such hardware is prohibitively expensive, which leads, perforce, to expedients, 
such as: (1) Leak testing with non-flammable tracer fluids (e.g. liquid Nitrogen); and (2) Leak testing with 
room temperature tracer fluids (e.y gaseous Helium). 

Such expedients undermine the realism of the tests. If, however, one could apply rational fluid dynamics 
methods to derive an general analytical expression with which one could relate the throughput of gaseous 
Helium through a given leak channel to the throughput of liquid Hydrogen through the same channel then 
one could recover much of the information that one would otherwise forfeit through these expedients. These 
facts lead to following questions: (1) What would be an example of a generic flaw in a gasket?; and (2) How 
can one calculate the flow of fluid it?. 

The following report addresses these questions. It considers a particular leak geometry, namely one 
formed by a gasket, a sealing surface, and a filament trapped between them (so that the cross section of the 
leak channel is a flat bottomed curvilinear triangle , two sides of which are circular arcs and which has cusps 
on all three corners). 
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LEAKAGE THROUGH A CHANNEL FORMED BY A GASKET, A SEALING 
SURFACE, AND A FILAMENT TRAPPED BETWEEN THEM 

John M. Russell 



Consider an apparatus of the form illustrated in Fig- 
ure 1, at left. One end of the metal bellows hose is 
plugged and the other end is attached to a residual gas 
analyzer capable of measuring small partial pressures 
of tracer fluid. The hose is evacuated and the plugged 
end is subject to any of a variety of test conditions, 
such as: (1) Immersion in a bath of liquid Nitrogen; 
(2) Exposure to cold gaseous Nitrogen by suspension 
above the liquid level in a liquid Nitrogen Dewar; (3) 
Exposure to gaseous Nitrogen at room temperature; 
(4) Exposure to gaseous Helium at room temperature. 
In this way, one may measure the leakage of tracer 
fluid that leaks from the outside of the plugged end of 
the hose to the inside. 

If the seal on the plug were flawless no tracer fluid 
would leak. One may, however, simulate the kind of 
flaw that may arise in practice by draping a fine wire 
or other hard filament across the nose seal prior to 
assembly (see Figure 1, where the filament is labeled 
“foreign object”). 

One would like to be able to model the flow 
through a leak channel of the above sort, but to do 
so, one must describe the shape of the resulting leak 
channel in detail. To this end, consider what a hypo- 
thetical observer might see if he or she could observe 
the filament along a line of sight parallel to the filament 
itself during assembly of the apparatus (see Figure 2 
on the next page). 

I have illustrated the spud and the plug as jaws of 
a vise. I have illustrated the nose seal in light shading. 
As the soft jaws close over the assembly, the soft seal 
material gives way to the hard filament as illustrated 
in Panel (b). One may surmise that the wire crevice 
will be effectively flat bottomed and that the arc of 


the crevice corresponding to the wire will be circular. There is no physical reason to assume that the arc 
corresponding to the gasket is also circular. The analysis simplifies considerably under this assumption, 
however, so I will make it in the remainder of this work. I will also use the term mire crevice to denote the 
kind of curvilinear triangle (with a flat bottom, two circular arcs, and cusps in the corners) illustrated in 


Panel (c) at the bottom of Figure 2. 

The remainder of this report is an analysis of the flow through a single wire crevice. Section 2 presents 
the necessary scientific background on the flow of liquids and gases through capillaries. Section 3 formulates 
the theory appropriate to the wire crevice. Section 4 presents the main result, namely a set of contours of 
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Figure 2. Generation of a wire crevice. Panel (a) represents the placement of the plug, nose seal, wire, 
and spud prior to assembly. Panel (b) represents these components after assembly. Panel (c) is a detail 
of the left member of the pair of leak channels illustrated in white in Panel (b). 

constant stre&mwise velocity (in nondimensional variables) for the flow in question, and Section 5 presents 
the conclusions and recommendations. 
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2. SCIENTIFIC BACKGROUND 

Authors of textbooks on fluid mechanics assign pride of place to the classical problem of fully-developed 
laminar flow through a pipe or channel and this problem applies to the present case when the fluid is a 
liquid. Many of the assumptions and conventions one introduces in the liquid-flow problem apply equally to 
the flow of a gas, so I will state them here with all due ceremony. 

Assumption of a prismatic boundary, orientation of coordinates and notation for velocity 
components and fluid parameters. One assumes that the boundary of the region occupied by the fluid 
has straight-line generators. Let (x, Y,z) be cartesian coordinates (each of which has the dimensions of 
length) with the positive z-axis parallel to the generators of the boundary surface and in the direction of the 
flow of fluid. Let (u, v, w) be the components of the fluid velocity vector (each of which has the dimensions 
of velocity) belonging to the coordinates (x, Y, z). Let P be the fluid pressure (a positive valued parameter 
with the dimensions of force-per-unit-area). Let p be the fluid density (a positive valued parameter with 
the dimensions of mass- per- unit- volume). Let p denote the dynamic viscosity of the fluid (a positive valued 
parameter with the dimensions of stress divided by strain rate, or, equivalently, force times time divided by 
area). 

Notation for typical lengths. Let L denote the diameter of the smallest circle into which the cross 
section of the channel can be fit. One may align the x-axis so that the two most widely separated points on 
a typical cross section are points on the X-axis and I will assume that the coordinate axes have this feature. 
Let 6 denote a typical dimension of the cross section in the direction of the Y-axis. The assumptions I have 
introduced thus far ensure that 6 < L. Finally, let t denote the streamwise length of the capillary. 

Slenderness assumption. In capillary flow theory, one assumes that L/£ = e 1. I will refer to the 
parameter e as the slenderness parameter in what follows. 

On the neglect of gravitational terms. One may interpret the differential equation for the rate of change 
of translational momentum of a fluid in the form of the statement “The acceleration of a fluid particle equals 
the resultant force-per-unit-mass (from all causes) to which that particle is subject.” One of the force-per- 
unit mass terms in the streamwise component of the momentum equation is the expression — (\j p){dp/dz), 
which represents direct driving by the unequal pressures at the inlet and the outlet. The presence of the 
minus sign reflects the fact that the fluid in a capillary flows in the direction of decreasing pressure. Another 
of the force-per-unit-mass term is — psin7, in which g > 0 is the acceleration due to gravity and 7 is the 
angle of inclination of the long axis of the capillary above the horizontal. Here 7 is positive or negative 
according to whether the fluid flows uphill or downhill, respectively. 

The neglect of the gravitational term, — #sin7, in comparison with the term due to direct driving by the 
pressure gradient, — (l/p)(dp/dz)> thus amounts to the assumption that the ratio of the orders of magnitude 
of these terms is very much less than one, viz. 


(l/p)\dp/dz\ <K 1 • (1) 

In the case of gas flow, we have the equation of state of an ideal gas, namely p = pRT, in which R is the 
gas constant particular to the gas in question (a parameter with the dimensions of energy per unit mass per 
unit temperature) and T is the absolute temperature. If one eliminates p from (1) by means of p = pRT , 
one obtains a restatement of (1) in the form 


p/\dp/dz\ 

RT/g 


( 2 ) 


One may estimate the numerator of the fraction in the left member of (2) by the capillary length £. 
The denominator RT/g has the dimensions of length and represents the change in elevation over which the 
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density of a stationary isothermal atmosphere decreases by a power of e. This so called scale height of the 
atmosphere is on the order of kilometers (about ten kilometers for air at 15° C). For a capillary whose length 
is on the order of millimeters, therefore, one may apply the assumption (2) with complete confidence. 

For a liquid, the assumption (1) amounts to the statement that the pressure difference due to the 
hydrostatic head of the working fluid between the inlet and the outlet is small compared to the applied 
pressure difference. If the applied pressure is on the order of atmospheres and the capillary length is on the 
order of millimeters, then this one may again apply (2) with confidence. 

Steadiness. I will restrict attention to the case when the velocity and pressure fields are independent of 
time. 

On the neglect of inertial reaction to streamwise acceleration and the exclusion of turbulent 
flow. A standard exercise in classical fluid mechanics is to estimate the ratio so called inertial to viscous 
terms in the momentum equation. If there were only one representative length scale A in the problem, then 
this ratio would be the Reynolds number Re = pWX/p , where W is a typical streamwise velocity. In 
capillary flow, however, there are two relevant length scales, namely the capillary length, £, and the smallest 
cross-stream length, 8 , of the channel. A measure of the inertial force- per- unit-mass is thus W 2 /i, while a 
measure of the viscous force per unit mass is p,W/(p6 2 ). One thus arrives at the following estimate of the 
ratio of the inertial to the viscous terms: 



The first factor in the right member is the slenderness factor, e, and the second factor is the Reynolds 
number, Re, based on the smallest cross-stream length of the channel. One may therefore restate the 
condition for the neglect of inertial terms in the momentum equation as follows: 


eRe 1 . (4) 

If i is on the order of millimeters and 6 is on the order of microns then e will be on the order of 10 -3 . In 
this case, Re need not be small and may even be in the hundreds without violating assumption (4). I will, 
however, confine attention to the case of non-tubulent flow, which generally requires that Re be less than 
2 , 000 . 

Assumption of constant temperature. I will assume that the temperature of the fluid (whether gaseous 
or liquid) is of uniform temperature. Since the fluid viscosity, p, is primarily a fuction of temperature, the 
constant temperature assumptions effectively guarantees that the viscosity is uniform throughout the fluid. 

The Poisson equation for the streamwise velocity. Every viscous term in the equation of motion of a 
fluid involves a second derivative of a velocity component with respect to some combination of the spatial 
coordinates. In view of the slenderness assumption e«l, one may separate the various viscous terms into 
classes distinguished by their respective orders of magnitude. The largest such terms are those in which both 
spatial derivatives are with respect to the cross-stream coordinates (x, y). Keeping these terms and making 
a few additional uses of the assumption one arrives at the following dominant balance of terms in the 

streamwise momentum equation: 

d 2 w d 2 w ^ 1 dp 
dx 2 + d\ 2 ~ /x dz ' 

in which P is a function of Z only. 

For a particular cross section z = constant, equation (5) takes the form of a POISSON equation for the 
distibution of streamwise velocity w across that station. The boundary condition for (5) is the condition of 
no-slip between a viscous fluid and an adjacent solid wall: 


w 


wall 


= 0 . 


( 6 ) 
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Introduction of nondimensional variables. One may introduce nondimensional cross-stream position 
coordinates (x, y) by the definitions 

x = x/L , y = y/ L . (7^,6) 

In the foregoing discussion, I have used the symbol W for a measure of the typical streamwise velocity across 
the cross section of the channel. One may make this definition more specific by taking 


WmX±. 

p dz 


(8) 


Note that fluid flows from greater to lesser pressure and the positive z-axis in the direction of the flow. It 
follows that dp/dz < 0, so the right member of (8) is positive as expected. Having a scale for the streamwise 
velocity, one may define the nondimensional velocity component w by 


w : 


w/W . 


( 10 ) 


If one multiplies the Poisson equation (5) by L 2 /W, one obtains, after some reductions, the following 
normalized form: 

+ ^ = (id 

3x ! V 


whose solution is subject to a boundary condition equivalent to (6), viz. 

= 0 . 


w\ 


wall 


( 12 ) 


The solution of the boundary- value defined by (1 1 ) and (12) will be a function ( x , y) ^ w. Note, in particular, 
that w does not depend upon anything other than x and y. 

Conservation of mass along the channel and its consequences. Let S' denote a surface of the x-v 
plane corresponding to the channel cross section. By the same token, let the unpunctuated symbol S denote 
the corresponding surface with respect to the nondimensional variables ( x,y ). One may write an expression 
for the rate of transport of mass, m, across a typical cross section of the channel by 


m = JJ pW dxdY . 

S' 


(13) 


Under the assumptions stated thus far, the mass density, p, is, at worst, a function of the streamwise 
coordinate z. One may therefore move the coefficient p outside of the integral sign in (13). If one carries 
out this step and also introduces the nondimensional coordinates {x, y) and the nondimensional streamwise 
velocity w by means of the definitions (7 a, b) and (10), one obtains 


m = pWL 2 



(14) 


If one eliminates W by means of (8) and rearranges slightly, one obtains 

rhfi dp 

L 4 JJ w dxdy ^ dz 
s 


(15) 
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The assumption of steadiness and the law of conservation of mass imply that m is independent of z. All 
of the other terms in the left member of (15) are also independent of Z under assumptions already stated. 
One concludes that for either liquids or gases 


9 


dp 

dz 


= independent of Z . 


(16) 


In the particular case of liquid flow, p is constant, so constancy of pdp/dz with respect to Z implies 
constancy of dp/dz . It follows that P varies linearly with respect to z and we have 


dp p 2 - pi 

P-TZ = P 

dz Z2 — z i 


P2 - Pi 
= P—f- 


(17) 


for liquid flow. Here, the subscripts “1” and “2” denotes conditions at the inlet to and the outlet from the 
capillary, respectively. 

In the particular case of gas flow, p varies according to the equation of state P — pRI\ so 


dp p dp 1 d(p 2 / 2) 

~ HTdz ~ 7 ?T dz 


Since the leftmost member of (18) is independent of Z (by (16)), one concludes from (18) and the assumption 
of constant temperature that p 2 varies linearly with respect to Z. It follows that 


dp _ 1 p| ~ P i _ l P 2 ~ Pi 

p dz 2RT Z 2 — z i 2 RT i 


for gas flow. 

One may combine the results expressed by (15), (17), and (19) in the following form 


f L J j w dxdy p^—j^—, 


m = < 


L 4 


^ Pi - p! 


,, jS wdxdv 2/rr t 


for liquid flow; 
, for gas flow. 


( 20 ) 


The two cases of equation (20) furnish an analytical correlation between the throughputs of liquid and gases 
through the same capillary. 

Authors of works on leak detection sometimes employ the symbol S to denote the rate of transport of 
fluid volume across a cross section of a fluid conduit. These same authors also sometimes write Q for the 
product pS and refer to the resulting quantity as the gas throughput. In view of the equation of state of an 
ideal gas, one may relate Q to rh by the following identities 


rh = pS = 



Q_ 

RT 


(21) 


If one multiplies (20) (in the case of gas flow) by RT , one obtains, in view of the outermost equality in (21), 
the identity 

( 22 ) 

s 
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The quantity Q so defined has the dimensions of pressure times volume transport rate. Thus, the expression 
atmcm 3 /s (also known as standard cubic centimeters per second) are suitable units for the measurement of 
Q. Of course, the outermost identity in (21) allows one to relate Q to the mass transport rate m. 

The case of a circular cross section; Classical formulas. The practical implementation of the foregoing 
formulas (which apply to all capillaries regardless of the cross sectional shape) depends upon the evaluation 
of the integral 

JJ w dxdy , (23) 

s 

which, in turn, depends upon the solution of the normalized boundary value problem defined by (11) and 
(12). The simplest example of a solution of this system is that of the capillary whose cross section is a 
circle of diameter L. In nondimensional coordinates, the diameter becomes one and, thus the radius becomes 
1/2. The equation of the boundary shape ( i.e . a circle of radius 1/2 centered on the origin) thus becomes 
(1/2) 2 - x 2 - y 2 — 0. If one notes that the the operator 


dx 2 ^ dy 2 


(24) 


takes the expression (1/2) 2 - x 2 - y 1 to a constant (in this case, -4) in the interior of the circle and that the 
expression (1/2) 2 —x 2 -y 2 vanishes on the boundary, it becomes plain that multiplication of (1/2) 2 - x 2 - y~ 
by a suitable constant (in this case, 1/4) yields a solution of the boundary value problem (11), (12) for the 
circular geometry. One concludes that 


(1/2) 2 - x 2 - y 2 


(25) 


in circular geometry. 

If one introduces (r, 6) for polar coordinates in the ( x,y ) plane, then (25) becomes 

(1/2) 2 - r 2 


( 20 ) 


and the integral (23) is equivalent to 




27 rr dr = 


7T 1 

8 ¥ ’ 


(27) 


If one substitutes this result into (20) and notes that p 2 — p 2 “ (pi +p2)(pi ~P2)> one obtains 


rn = 


JL ( iiV Pi ~P2 

< 8/*W P *■ ’ 

JL ( L\ 4 ( (Pl+P 2 )/ 2 ^ Pi ~P2 

\ RT ) e 


for liquid flow; 
for gas flow. 


(28) 


The liquid-flow form of (28) is associated with the names of Hagen and Poiseuille, both of whom studied 
the flow experimentally in the nineteenth century. The result was derived analytically by George Gabriel 
Stokes. The earliest analytical derivation of the gas-flow form of (28) appears to be that of Oskar Emil 
Meyer in 1866 (Ref. 1, page 269). 
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In view of the ideal gas law, the constant valued expression 

(p i+P2)/2 
FT 

has the dimensions of density, and one may regard it as an effective constant density for the gas in a 
capillary (based on the ideal gas law and the arithemetic mean of the inlet and outlet pressures). While 
this interpretation is convenient as a mnemonic device, it is important to realize that the derivation of the 
Meyer formula does not assume constancy of the density even as an approximation. Indeed, the Meyer 
formula holds for any expansion ratio. The foregoing derivation supports the last assertion nothwithstanding 
contrary statements by modem authors in the literature on leak detection. 


3. FLOW THROUGH A WIRE CREVICE 

IVansformation to a Dirichlet problem. Let (x, y) w p be a function defined by the equation 


u>p = -y 2 / 2 . 


(29) 


It follows that 


d 2 w p d 2 w p 
dx 2 dy 2 


(30) 


One concludes that w p is a particular solution of the differential equation (11) satisfied by w. Let (x, y) >—> <j> 
be a function defined by the equation 

i u=w p + <t>. (31) 


If one subtracts (30) from (11), and uses the identity w -w p = <f> (which follows from (31)), one obtains 


d 2 <fi d 2 <j> 
dx 2 ^ dy 2 


(32) 


If, now one substitites (31) into the boundary condition (12), one obtains 

wall “ I wall “ I wall 1 


or 



(33) 


in which the subscript ( ) w denotes evaluation on the wall. 

The partial differential equation (32) that <j> satisfies is Laplace’s equation. A function that satisfies 
Laplace’s equation is called a harmonic function. The boundary- value problem consisting of Laplace’s 
equation and a statement of the boundary values of the unknown (as distinguished from the normal derivative 
of the unknown or some combination of the unknown and its normal derivative) is called the Dirichlet 
problem. The pair of equations (32) and (33) thus constitutes a DIRICHLET problem for <f>. 


Construction of harmonic functions from differentiable functions of a complex variable. Expe- 
rience has shown that one may expedite the solution of the DIRICHLET problem in two space dimensions 
considerably through the use of complex variables. There are two reasons why this is so. To illustrate the 
first reason, let £ = f + it) be a complex number. Here t = y /^ T and the pair of real numbers (f , r?) represent 
a point in a plane (the complex C-pl&ne). If C F is an analytic ( i.e . differentiable) function that maps C 
to the complex number F (with the expansion, say, F = <f> T then the real and imaginary parts of F 
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(namely <p and xp, respectively) are related to the real and imaginary parts of £ (namely £ and r/, respectively) 
by the Cauchy-Riemann equations 


d<p dip d<p dip 

d£ dr] ’ dr] d£ 

If one eliminates xp from the (34 a, 6) by cross differentiation, one obtains, immediately, 

at? a n 2 


(34 a, it) 


(35) 


In this way, any analytic function of a complex variable ( e.g . C F) is related in a systematic way to a 
harmonic function (e.g. (f, r?) <p). 

A class of changes of independent variable that leaves Laplace’s equation invariant. A second 
reason why complex varibles are useful in solving the Dirichlet problem in two dimension concerns the 
invariance of Laplace’s equation under a class of transformations of the independent variables. Suppose, 
for example, that z »-* m is an analytic function that takes the complex number 2 to the complex number 
m and that the expansions of 2 and m into real and imaginary parts are 


z = x -f iy , rn = m r T irrii . 


(36) 


Then the function z i-+ m induces a transformation (x,y) (m r ,mi) (called a conformal transformation) 
between points in the complex 2 -plane and points in the complex m-plane. A basic exercise in the theory of 
functions of a complex variable is the demonstration that 


d 2 d 2 

dx 2 + dy 2 


dm 2 / d 2 d 2 \ 
dz \ dm 2 ^ dm 2 ) 


(37) 


(cf. Jeffreys & Jeffreys, p534) (Ref. 2). If, for example, (x, y) <f> is harmonic, then 


d 2 <p d 2 (j> 

dx 2 + dy 2 


dm 2 / d 2 cj) d 2 <p \ 

dz \ dm 2 dm 2 ) 


(38 J 


Division of the outermost equality by \dm/dz \ 2 yields 


d 2 <f> d 2 <j> 

dm 2 + dm 2 


(39) 


Thus, if </> is a harmonic function of one set of independent variables, and those variables are related to a 
new set of variables by a conformal transformation then <}> will also be harmonic with respect to the new 
variables. 

TVansformation of the wire crevice to the interior of the unit circle. The book Analytic function 
theory by Einar Hille (Ref. 3) gives several examples of conformal transformations, of which two are 
particularly germane to the present problem. Figure 3 illustrates regions in the planes of four distinct 
complex variables. The transformation between what I have called the m-plane and the r-plane appears on 
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Figure 3. Three step conformal transformation of a wire crevice (in the 2 -plane) to the interior of a unit 
circle (in the £-plane). 


(and near) page 177 of Hille’s book and the transformation between what 1 have called the r-plane and 
the £ plane appear on (and near) page 370 of that source. 

HiLLE’s formula for the transformation m r is equivalent to 


itf(l-m) 

K(m) 


(40) 


in which m K(m) is a so-called complete elliptic integral of the first kind 
is 



(ibid., pl50). The inverse of (40) 

(41) 
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in which (z,r) »— ► di(z\r) ( i 6 {1,2, 3,4}) are the so-called theta functions , (ibid., pl56). Some modern soft- 
ware packages, such as Mathematica , enable evaluation ofm^ K(m) and (z y r) ■— ► t9i(^|r) (i € {1,2, 3,4}) 
for complex arguments as easily as the evaluation of sines and cosines. The use of such higher functions 
is therefore not a practical limitation on a the development of a computational algorithm for the present 
problem. 

Hille’s formula for the transformation ( m m (or a trivial modification therof) is equivalent to 


•m- 


(42) 


whose inverse is 


< 


= i( — ) ■ 

\m + i ) 


(43) 


HlLLE (ibid., p 46, et seq.) discusses a class of conformal transformations which he associates with the 
name Mobius. In these transformations, two complex variables are related by a formula which equates the 
first variable to a ratio of linear functions of the second (The inverse of such a transformation equates the 
second variable to a ratio of linear functions of the first). HlLLE calls attention (ibid., page 51) to the class 
of geometric figures consisting of circles and straight lines. A basic feature of the Mob I US transformation is 
that it takes a circle or line in the plane of the first complex variable to a circle or line in the plane of the 
second. In this sense, the class of circles and lines is invariant under Mobius transformations. With this 
theorem in mind, it is not difficult to find (by trial and error) a MOBIUS transformation that takes the wire 
crevice in the z-plane to the interior of the shaded figure in the r-plane. Such a Mobius transformation is 
given by 

1 Z j Z CKlsp 


whose inverse is 


1 - 2 


r — 1 


T {\/ z cusp) 


(44) 


(45) 


Here z cusp is the intersection of the gasket and the wire in the wire crevice problem. 

The condition that the gasket and the wire be tangent to each other restricts values of z cusp to points 
on a semicircle of diameter one centered at the point 1/2 on the real axis. One may write a formula which 
encompasses all possible values of z cusp that fulfill this condition in the form 


*xusp 


1 

l-iA 


(46) 


in which A > 0 is a real parameter. In fact A 2 = R^/R gy where and R g are the radii of curvature of the 
wire arc and the gasket arc, respectively. 

Boundary value of <f> in the £-plane. The foregoing conformal transformation formulas enable one to 
associate any point ( x , y) in the in or on the wire crevice in the 2 -plane (where z = x + iy) to one and only 
one point (f, rj) in or on the unit disk in the C-plane (where C = f + «?)• Thus, the boundary values (33) for 
the unknown harmonic function (x, y) •— * <f> in (32) carry over to corresponding boundary values of in the 
C-plane. Figure 4 illustrates the corresponding distribution of these boundary values expressed in terms of 
the polar coordinate angle 6 in the C-plane (The reader should note that my use of (r, 6) as polar coordinates 
in the C-plane here is a new use of these symbols, not to be confused with my earlier use of them in the text 
between (25) and (26) above). 
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4>w 



Figure 4. Distribution of the boundary value of the harmonic function <p as a function of polar angle in 
the plane 


Solution for <f> in the £-plane as an infinite series. One may relate the Laplace operator in cartesian 
coordinates (£, 77 ) in the £-plane to polar coordinates (r, 6) in that plane by the identity 


0 2 ( ) , a 2 ( ) 

3£ 2 dr) 2 


ld_ 

r dr 


d()\ , 1 d\) 


m 


+ 


dQ 2 


(47) 


One may verify that for any integer n, the operator in the right member of (47) takes expression r n cos(n0) 
to zero. The same operator takes r n sin(ntf) to zero and one arrives at the conclusion that any function of 
the form 

00 

4> = 7 T + ^2 r ”Ia„ cos (n0) + b n sin(n0)] , (48) 

n— 1 

whose constant, real- valued coefficients (on,fc n ) are subject to the condition that the series (48) converge, is 
a solution of Laplace’s equation in the C-plane. In view of the foregoing discussion under the heading “A 
class of changes of independent variable that leaves Laplace’s equation invariant” such a representation for 
<p must also satisfy Laplace’s equation in the z-pl&ne. 

In view of the complex variable identities cos0 — ( e t& + e~ t9 )/2 and sin 0 = (e t0 — e~ l$ )/(2 1 ), one may 
arrange (48) in the following equivalent form 


4> = 





(49) 


One may shorten the notation by defining values of a n and b n for zero and negative subscript as follows ( e.g . 
LlGHTHILL (1959, p3), Ref 4): 

U— n = O n , 6 _ n = b n (50) 
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(so &o = 0, by construction). Moreover, for any integer n, we write 

(51) 

With these conventions, equation (49) takes the form 


oo 

<j> = J2 Cnr M e in9 . 

n~~ oo 


(52) 


If one substitutes r = 1, one obtains the boundary value of <f>, namley 4> w , so 

CO 

4> w = Y, c " eine ■ ■ {5:}) 

n= — oo 


I have already remarked that the boundary- value function 0 *-* <f> w is known (c/. Figure 4 above). Let k be 
an arbitrary integer. If one multiplies (53) by e~ ikd and integrates with respect to 0 from zero to 2n y one 
obtains 

r 2n oo r 2 n 

/ <t> w e- ike d9 = Y' C n / e i{ - n ~ k)e dO . (54) 

Jo Jo 


But 


It follows that 


t (n-k)0 __ 


r2n 

/ e i{n - k)e dJ9 

Jo 

If one substitutes this last result into (54), one obtains 


JL( 

dO ' 


>) 

if k - n\ 

e i(n-k)0' 

j if k j-- n. 

i(k—n) 

f 2t r 

if k — n; 

= \0 

if k n. 


(55) 


(56) 



<j> w e- ik0 dO = 27rC fc , 


(57) 


or, equivalently, 



(58) 


(which holds for all integer n. Equation (58) thus determines the coefficients in the series representation 
(52). 

Solution for w in the z-plane. Knowing 4> at any point in the C-plane, one may associate the same 
value of 4> with the corresponding point in the z-plane ( via the conformal transformations given in equations 
(40)-(45)). One concludes that </> is now known through the whole wire crevice. Since w p is specified by 
(29), one concludes that the whole right member of (31) for tv is known throughout the wire crevice. 


4. RESULTS 

Figure 5 gives a plot of the contours of constant w obtained from a numerical evaluation. The slight 
raggedness of the outermost contour is probably due to an overly course discretization of 8 in the calculation. 
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Figure 5. Contours of constant streamwise velocity in the flow through a wire crevice. The innermost 
contour represents 3/4 of the maximum value of w, the next innermost contour represent 1/2 the 
maximum, and so on. 

5. CONCLUSIONS AND RECOMMENDATION 

The foregoing work supports the following conclusions: 

1. The liquid-gas flowrate correlation embodied by (20) (which holds for arbitrary shaped cross sections) 
and (28) (which holds for circular cross sections) are on a firm scientific basis and are suitable for 
immediate application. 

2. One may solve the boundary- value problem for the flow through a wire crevice analytically and the 
solution is expressible in terms of a single infinite series. 

3. The solution for wire-crevice geometry is computationally feasable on a small computer. 

I had hoped to refine the calculations so as to eliminate the raggedness of the tc-contours before com- 
pletion of the summer term, but ran into some snags during the last two weeks. The recommendation, of 
course, is to continue the effort including the calculation of the integral (23), and thence, the throughput. 
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